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Exact theoretical results for the violation of time dependent fluctuation-dissipation relations in 
driven dissipative systems are presented. The ratio of correlation to delayed response in the stochas- 
tic model introduced in [Phys. Rev. Lett. 93, 240601 (2004)] is shown to depend on measurement 
time. The fluctuation temperature defined by this ratio differs both from the temperature of the 
environment performing the driving, and from other effective temperatures of the system, such 
as the average energy (or "granular temperature"). General explanations are given for the time 
independence of fluctuation temperature for simple measurements or long measurement times. 
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The fluctuation-dissipation relation (FDR) provides a 
fundamental connection between two small deviations 
from thermodynamic equilibrium - the autocorrelation 
function of an observable, and the response function of 
that observable to changes in its conjugate field. In equi- 
librium, there are two aspects to the FDR: (i) these 
functions have the same spatial and temporal (or fre- 
quency) dependence, and (ii) for every observable their 
ratio equals the system's temperature. 

Away from equilibrium, there is no such theorem con- 
necting correlation and response. This having been said, 
there is much theoretical 0,|llilllllilli,|i|lS|ll|l2 
and experimental |l.'lll 1| interest as to what the relation 
between these quantities is. In particular, one may ask 
whether either aspect (i) or (ii) holds for small deviations 
from non-equilibrium steady states (NESS). Aspect (ii) 
has been tested extensively for low frequency measure- 
ments in various systems, and the effective temperature 
Tp, defined by the ratio of correlation to response, clearly 
differs from the environment temperature [ij, |2j, Il3| . On 
the other hand, numerical simulations indicate that Tp 
does coincide |15| in some cases with other effective tem- 
peratures defined for glasses , sheared dissipative foam 

, dense granular packings |J and driven granular gases 

Regarding aspect (i), different frequency dependence of 
the correlation and response functions has been observed 
in simulations of glasses Q , sheared foam Q and models 
of driven systems , as well as in bacterial bath exper- 
iments jjjj • In ah these cases, fast modes thermalizc, 
while slower modes manifest a higher (usually frequency- 
independent) temperature reflecting driving or system 
history (e.g., temperature quench). Experiments on vi- 
brated granular systems |14| have shown a weak depen- 
dence of Tp on measurement frequency, which may not 
obviously be attributed to thermalization of fast modes 
with the environment. 

This work aims at understanding the relation be- 
tween time- (or frequency-) dependent correlations and 
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response in NESS. To our knowledge, these have not pre- 
viously been calculated exactly for any driven dissipative 
system. Here we calculate the autocorrelation and re- 
sponse functions in the context of a simple stochastic 
model of a dissipative system. The model introduced 
in [l2j is constructed about the essential features of any 
driven dissipative system: energy loss upon interaction 
between particles (or modes) and energy input through 
an external driving mechanism. As such, we have reason 
to hope that our results will be relevant to a broad class 
of driven dissipative systems . 

Our main results are: (1) There are certain observables 
for which correlation and response have identical tempo- 
ral dependence, however, (2) For a general measurement, 
they have different temporal dependences. This notwith- 
standing, (3) For asymptotically long times, the time de- 
pendence of correlation and response coincides. Spatial 
FDR in our model are investigated in [18j . 

The model consists of N particles, each with two de- 
grees of freedom (DOF), one "kinetic", e.j, and one "in- 
ternal", Xi. The energy of particle i is e, — X{Fi, where 
Fi is an external field on particle i, which may depend 
on time. The {xi\ may be thought of as positions or any 
other coordinates used to measure FDR. In our notation 
we emphasize the difference between microscopic DOF's 
{ei} and {xi} and macroscopic quantities E = Yli=i e * 
and X = X)j=i x i- The system dissipates energy through 
interactions and is maintained in a steady state by being 
coupled to a thermal bath. 

In every interaction two DOF's are randomly chosen 
and their energy is stochastically redistributed between 
them. In interactions between two "kinetic" DOF's 
(e — e) there is dissipation: only a fraction < a < 1 
(a is like a restitution coefficient) of the energy is con- 
served and the remaining is dissipated out of the system. 
Interactions between a "kinetic" and an "internal" DOF 
(e — x), and between a "kinetic" DOF and a DOF from 
the bath (e — eg) conserve energy. When all interactions 
are elastic (a — 1), detailed balance holds and the system 
reaches thermodynamic equilibrium with the bath. Rates 
of interactions are determined by the per-particle inter- 
action rate T, and the dimensionless coupling strengths 
< (/, h) < 1 of the "kinetic" DOF's to the bath and to 
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the "internal" DOF's, respectively. 

The resulting stochastic equation of motion for particle 
i reads, 



value : 



probability : 



e t (t + dt) = 



ei(t) 1 - Tdt 

za[et(t) + e k (t)] (1 - f)Tdt (la) 

z[ei(t) + e B ] f(l-h)Tdt 
z[ei(t) - Xji^F^t)] fhTdt 



value 



Xi(t + dt) 



Xi(t) 

- ei (t)/fi(t)] 



probability : 
1 - fhTdt 
fhTdt 



(lb) 



where: j, k G {1, N} are indices of particles with which 
particle i may interact, chosen randomly at every interac- 
tion (we do not allow e^ — interactions, therefore k ^ i); 
z e [0, 1] is the fraction of repartitioned energy given to 
particle i in the interaction, chosen randomly from a uni- 
form distribution at every interaction; es is the energy of 
a bath DOF chosen randomly at every interaction from 
the equilibrium distribution pb(^b) — Tg 1 e~ eB / TB at 
the environment temperature Tg. 

In the steady state we define 1-particle and N-particle 
autocorrelation functions, c(t) = (xi(t)xi(0)} — (xi) 2 and 
C{t) = (X(t)X(0)) — (X) 2 , respectively, and correspond- 
ing delayed response functions r(t) = d(xi(t)) /dFi(0) 
and R(t) = d(X(t))/dF(0). The responses r(t) and R[t) 
are to a sudden change at time t = of the field on the 
single particle i or of the uniform field on all particles, re- 
spectively. ( ) denotes the steady state ensemble average 
over possible states at the beginning of the measurement 
(t = 0). In equilibrium (a = 1 or / = 1), c(t) = r(t)Ts 
for any x$, and C(t) = R(t)Ts for X |l9j; we will be 
concerned with the non-equilibrium cases. 

For 1-particle FDR we consider a system evolving from 
a given state at t = until a later time t. We average 
Eq. HQ over the stochasticity in the dynamics for this 
initial state, and denote this averaging by an ovcrlinc, 



2 dxi (t) 

r dt 



-fh Xi (t) 



E(t) 
NF{t) 



(2) 



then, at time t = 0, changing the field Fi actin g on 
particle i from Fq to F. We then follow Xiit) and 
average over the system's steady state distribution at 
time t — to yield (xi(t)). The steady state solution 
of Eq. © yields -(x t (0))F o = -{x t )F = (e), thus 
(xi(t)) = — (e) [F^ 1 + (Fq 1 - F-^e-i*]. After differen- 
tiating with respect to Fq and taking the limit Fq — * F, 
we get r{t) = (e)F- 2 e-"< t . 

We thus see that the 1-particle correlation and re- 
sponse have the same temporal dependence and aspect 
(i) of the FDR exactly holds, with the time- independent 
fluctuation temperature plj . 



rpl _ C(t) 



Tr 
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Aspect (ii), on the other hand, is violated [T^ , since Tp 
generally differs both from the environment temperature 
Tb and from the granular temperature, defined as the 
average energy per DOF: Tq = (e). Only in the equilib- 
rium limits (a — > 1, / — > 1) is the distribution of X\ the 
exponential Boltzmann distribution and Tp = Tq = Tb- 

The FDR is only proven for equilibrium, so one might 
expect the correlation and response to depend differently 
on time in a NESS. However, the correlation represents 
the system's return to steady state after deviating from 
it due to spontaneous fluctuations, and the response rep- 
resents its return to steady state after being moved away 
from it by some external force. Both are governed by the 
same physical processes and thus generally posses the 
same time scales. In simple cases (as the 1-particle mea- 
surement solved above) there is only one time scale and 
the correlation and response are not rich enough to have 
different time dependences. In order to observe violations 
of aspect (i) in the FDR we now consider a measurement 
on a macroscopic quantity. This possesses two time scales 
- just enough to exhibit a different temporal dependence 
for the correlation and response. 

We assume uniform F and sum Eq. J2J over particles, 



2 dX(t) 
f dt 



= -fh X(t) 



E{fy 
F(t) 



(5) 



By the central limit theorem, for large systems (N 3> 
1) the relative fluctuation (xi — (xi))/(xi) of xi is much 
larger than the relative fluctuation (E — (E))/{E) of E, 
thus the steady state average (E) may be substituted for 
the instantaneous value E(t) in Eq. (0). For a constant 
field, notingthat (e) = (E)/N = -\xi)Fi this yields, 



;(*) = (Xi) + (Xi(0) - (Xi))e 



-7* 



(3) 



with 7 = fhT/2. The correlation is obtained by multi- 
plying Eq. by Xi(0) and averaging over the steady 
state distribution for all possible initial (t = 0) states of 
the system, which yields c(i) = {(xf) — (a;i) 2 )e _7t |2lj . 

The 1-particle response is obtained by maintaining 
the system in a steady state with some external field, 



Now, fluctuations in E are not negligible compared to 
fluctuations in X, so Eq. JSJ must be solved in conjunc- 
tion with the equation similarly derived from Eq. Hla(l . 



2 dE(t) 

r dt 



-fhX(t)F(t) - A x E{t) + /(l - h)NT B , (6) 



where A n = n + (1 — /)(1 — 2ct n ). The average evolution 
of the system from a state with X(0) and E(Q) is hence 
given by the simultaneous solution of Eqs. JSJ and ©: 



X(t) = (X) 

2 r 



E 



a e (X(0)-(X)) + b t 



E(0) - (E) 



F 



(7) 
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1 /2 

where we denote k = (A\ - 2A x fh + 5f 2 h 2 ) ' , at = 
(1 - A 1 + fh)/2k, a 2 = (l + A 1 - fh)/2k, b x = fh/k, 
b 2 = -fh/k, 7i = (A 1 + fh + fc)r/4 and 72 = (A x + 
fh — fc)r/4, and the steady state values are given by 

-{X)F = (E) = /(l - h)NT B /(A 1 - fh). 

To compute the correlation we multiply X(t) by X(0) 
and average over the steady state distribution of the ini- 
tial states (that is, over X(0) and E(0)) 21]: 



at) = J2h((x 2 )-(x)' 

. (XE)-{X){E) 



F 



-"tit 



(8) 



For the response we change the field from Fo to F 
at t = 0, and after averaging Eq. Q) over the initial 
states taken from the steady state corresponding to the 
field Fo, differentiating with respect to Fo and taking the 
limit Fq — > F we have, 



(9) 



In this N-particle measurement C(t) and R(t) share the 
rates 71 and 72, however with different prefactors, leading 
to violation of aspect (i) in the FDR [22] • If one insists on 
defining a fluctuation temperature Tp (t) = C{t)/R{t), it 
will depend on the measurement time, as can be seen in 
Fig. ^ which also demonstrates that the 1-particle and 
N-particle measurements yield fluctuation temperatures 
differing both one from the other as well as from the 
granular temperature Tq and bath temperature Tb- 

As the correlation and response generally share the 
same time scales, for measurement times longer than the 
maximal time scale of the system (I/72 for the N-particle 
measurement in our model) this time scale dominates 
both correlation and response and aspect (i) of the FDR 
asymptotically holds for long times (see Fig. It is 
interesting to speculate that this behavior may hold for 
general dissipative systems. 

Although the FDR is generally not valid in NESS, it 
is hard to observe its violations in simulations or experi- 
ments on dissipative systems. One reason for this is that 
dissipative systems often have nearly Boltzmann distri- 
butions with some effective temperature, and deviations 
between different effective temperatures of the system 
(like Tp and Tq) are too small to be observed 0,03 • 

Another reason is that some measurements on some 
systems do not exhibit at all the FDR violations pre- 
sented above. To demonstrate this we now show that 
in the inelastic Maxwell model the frequency depen- 
dent Kubo formula 24] holds exactly with an effective 
temperature coinciding with Tq. That is, the velocity 
autocorrelation D(u) = f °° (v(0)v(t))e~ lujt dt relates to 
the mobility n(ui) by D(uj)= /j,(u)Tk with Tk = Tq in- 
dependent of frequency [2£j . We shall first show that for 
our model Tk is frequency- independent and equals Tp, 




R(t)/N 
C(t)/N 



C(t)/R(t) 



c(t)/r(t) 



FIG. 1: (color online) (a) Correlation and response functions 
and (b) resulting fluctuation temperatures vs. measurement 
time for the model with a — 0.6, / = 0.7 and h = 0.9. Graphs 
are scaled to Tb = 1, F = 1 and V = f . The resulting decay 
rates are 7 = 0.32 , 71 = 0.72 and 72 = 0.07. Numerical sim- 
ulations with = 10 3 (dashed lines) agree with theoretical 
results for N ^> 1 (thick gray lines). The bath temperature 
Tb and the granular temperature Tq are given for reference. 



which generally differs from Tq, and subsequently show 
that for the Maxwell model Tk — Tq. 

To make contact with the standard Kubo relations 
we will imagine the DOF's {xi\ in our model as po- 
sitions; thus the autocorrelation of the single particle 
"velocity" vi = dxi/dt is obtained by twice differenti- 
ating c(t) calculated above. The relation c(t) = r(t)Tp 
with Tp independent of t, as obtained for the f-particle 
FDR solved above, may be twice differentiated with re- 
spect to time and transformed to frequency domain. The 
Kubo formula then immediately follows with Tk = Tp. 
We have also verified this general result by explicitly 
calculating D(uj) and /x(w) in our model and obtained 
D(w) = n{u)T F = - (x,) 2 )/( 7 + iw). 

Our method may now be applied to the ID inelastic 
Maxwell model (the 2D and 3D versions follow trivially) : 
During a finite time dt > 0, short compared to the colli- 
sion rate T and to the coupling rate A to the thermostat, 
the effect of the fluctuating force in the Langevin dy- 
namics each particle undergoes is proportional to Vdt 
(see e.g. [2^). Thus the velocity of particle i evolves as, 



val 



ue 



probability 



Vi(t + dt) 



(1 - Xdt)vi(t) +ipi{t)Vdt 1-Tdt (10) 
^Vi(t) + i^v/(t) Tdt 



For stochastic thermostating, ^i(t) is an uncorrelated 
random force with (ip 2 ) = 2XTb, while for Gaussian ther- 
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mostating, ipi(t) = and A = (a 2 — l)T/4. Averaging Eq. 
(PI) yields (v(O)v(t)} = (v 2 )e~ Kt with k = A+(l + a)r/2, 
thus D{uj) = (v 2 ) /(n+iui). In the mobility measurement, 
a periodic acceleration £(t) — £oe Mt is added to a single 
particle, and from Eq. lfTU|) we have 



^ = - K ^ + &e -*. (11) 

This has the steady solution (uj(t)) = £oe tuJt /(« + «*>), 
thus //(w) = l/(/s + iw), and the Kubo formula holds 
with Tk = Tg- 

This difference between our model and the Maxwell 
model derives from the fact that Tp and Tk measure 
fluctuations, expressed in second moments of DOF's. In 
the Maxwell model Tq is defined from (v 2 ), which is the 
second moment of a DOF, while in our model Tq = (e) 
is defined from the first moment of a DOF. Only in equi- 
librium do all moments of the exponential Boltzmann 
distribution yield a single energy scale and these mea- 
surements in both models give the actual temperature. 

To our knowledge, we have presented the first exact 
results for time dependent violations of FDR in driven 



dissipative systems by demonstrating that correlation 
and response functions have different temporal behav- 
iors. This violation, like differences between different ef- 
fective temperatures of a system, may be small, requiring 
sensitivity in simulations and experiments. We identified 
situations relevant to general driven dissipative systems 
where the correlation and response have the same time 
dependence: In sufficiently simple measurements the cor- 
relation and response share a single time scale and so have 
the same dependence on time, thus time-dependent FDR 
violations may be observed only in systems or measure- 
ments exhibiting multiple time scales. Finally, since the 
correlation and response typically share the same time 
scales, their temporal behavior for measurements with 
long waiting times asymptotically coincides. 
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